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The link between the inversion gain of quantum cascade structures and the Bloch gain in periodic 
superlattices is presented. The proposed theoretical model based on the density matrix formalism 
is able to treat the gain mechanism of the Bloch oscillator and Quantum cascade laser on the same 
footing by taking into account in-plane momentum relaxation. The model predicts a dispersive con- 
tribution in addition to the (usual) population-inversion-dependent intersubband gain in quantum 
cascade structures and - in the absence of inversion - provides the quantum mechanical description 
for the dispersive gain in superlattices. It corroborates the predictions of the semi-classical miniband 
picture, according to which gain is predicted for photon energies lower than the Bloch oscillation 
frequency, whereas net absorption is expected at higher photon energies, as a description which is 
valid in the high-temperature limit. A red-shift of the amplified emission with respect to the reso- 
nant transition energy results from the dispersive gain contribution in any intersubband transition, 
for which the population inversion is small. 



I. INTRODUCTION 

Soon aftec. the original proposal of semiconductor 
superlatticesEJ, two apparently quite different schemes to 
obtain optical gain in such novel systems were put for- 
ward. In a seminal work, Kazarinov and Suris pointed 
out how to achieve a population inversion, a key ingredi- 
ent to obtain light amplification, between electronic sub- 
bands in a strongly biased superlatticco. On the other 
hand, based on semi-classical arguments, Ktitorov et alu 
and later Ignatov et al.U predicted optical gain due to 
Bloch oscillations within a miniband - despite a missing 
population inversion. 

Two decacks later, the demonstration of the quantum 
cascade laserB affirmed the first proposal. The (conduc- 
tion) band structure in each period is carefully designed 
to allow for injecting electrons into an upper subband 
state, with a long non-radiative lifetime, and to enable 
a fast extraction of electrons from an accordingly tai- 
lored lower state. As a consequence, population inver- 
sion is achieved. Also, by a suitable design, the struc- 
ture is electrically stable at threshold. By now, the 
quantum cascade laser technology covers a wide range 
of the electro-magnetic spectrum.— Recently, a room- 
temperature continuous-wave laseifl emitting at a wave- 
length of 9 fim has been demonstrated and stimulated 
emission in the terahertz regime at about 66 /imLnj has 
been observed. 

In contrast to this, the feasibility of the Bloch oscilla- 
tor, emitting electro-magnetic radiation, tunable by the 
external electric dc field, is still under question. Besides 
the task to stabilize the electric field domains in a biased 
superlattice at the point of operation, the description of 
the gain mechanism is, so far, based on semi-classical 
models only. In a naive picture, electro-magnetic radi- 
ation of photon energy hu> w Tiujb = eFd is expected, 
corresponding to the frequency of Bloch oscillations uib 
that linearly depends on the applied dc field F. For, e.g. 
a superlattice period d of some nanometers and fields of 



several tenth of kV/cm, the photon frequencies are in the 
terahertz range. 

In fact, semi-classical calculations exhibit neither gain 
nor absorption at resonance, i.e. for to = u>b, but trans- 
parency. Particularly in the quantum-mechanical pic- 
ture, it is evident that only spontaneous transitions can 
occur at resonance. At sufficiently high electric fields, the 
miniband is split into the Wannier-Stark ladder, a set of 
states evenly spaced by eFd in energy. Resonant stimu- 
lated emission processes between adjacent states are bal- 
anced by absorption processes, because of the transla- 
tional symmetry of the system, which dictates equal oc- 
cupation for all rungs of the Wannier-Stark ladder. 

Nevertheless, the semi-classical calculation does pre- 
dict gain - without inversion - for non-resonant transi- 
tions with a "too small" photon energy, hu < fiujb, and 
absorption for tuv > hub- But the existence and strength 
of this semi-classically predicted Bloch gain is still dis- 
cussed, despite experimental observation of many re- 
lated phenomena in-superlattices, such as negative differ 
ential conductivity^, the associated Bloch oscillationsli2l 
and the .coMeling of the superlattice to external THz 
radiationO'EZl, just to mention a few. In particular, the 
gain mechanism is lacking an interpretation in terms of 
the quantum-mechanical Wannier-Stark picture. 

In this paper such a quantum-mechanical interpreta- 
tion of the Bloch gain in superlattices is suggested and 
a link is established between the intersubband gain orig- 
inating from a population inversion, with its symmet- 
ric spectral shape centered at the transition energy, and 
the dispersive gain predicted for a periodic superlattice, 
with its nearly anti-symmetric profile. The quantum- 
mechanical model, based on the density matrix formalism 
similar to the one employed earlierli3, yields a general ex- 
pression for the gain profile in intersubband transitions. 

In Sec. II we present the model system, and discuss as- 
sumptions and details of the density matrix calculation. 
With less stringent approximations than those made by 
Kazarinov and Suris we find an expression for the coher- 
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ence between two, at first spatially separated, subband 
states that are coupled by tunneling and broadened by 
intra-subband scattering. The coherence determines cur- 
rent density as well as optical transitions. Transforming 
to the basis of eigenstates of the biased heterostructure - 
the Wannier-Stark basis for the supcrlattice - the model 
is extended to describe optical transitions between any 
pair of subbands. 

In Sec. Ill we apply the theory to superlattices and ob- 
tain the quantum mechanical counterparts to the semi- 
classical results for both the Esaki-Tsu current-voltage 
characteristics and the dispersive Bloch gain. The re- 
sults are quantitatively compared to the predictions of 
the semi-classical picture, where good agreement is found 
for higher electron temperatures. 

In Sec. IV intersubband transitions are investigated for 
the quantum cascade structure and the relation between 
(anti-symmetric) Bloch gain and Lorentzian-shaped in- 
tersubband gain becomes apparent: the theory predicts 
a transition from the Lorentzian shaped inversion gain to 
the dispersive Bloch gain with decreasing population in- 
version accompanied by a red-shift of the peak gain with 
respect to the transition energy. 



II. THEORY 

We consider two subbands, confined in adjacent wells, 
that serve as a model system for photon-assisted tunnel- 
ing structures and, in particular, for transitions within 
the Wannier-Stark ladder of a weakly coupled superlat- 
tice. To start with, the same basis as in the original work 
of Kazarinov and Suris is chosen as an appropriate basis 
set. The wave functions \ik) are given by the product 
of the envelope functions ^i(z), maximally localizedlil in 
well i, and plane waves 



(r,z\ik) = %(z)e 



ik-r 



(i) 



where the z-axis is defined by the direction of growth, 
k = {k x ,k y ) = k denotes the in-plane momentum and 
r = (x, y) the lateral position. The matrix elements of 
the Hamiltonian in this basis are given by 



K k , = {ik\n\jk') = H%s kk 



y kk' ' 
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where the respective contributions H and V take the 
form 
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dispersion relation parallel to the layers in the effective 
mass approximation 



n 2 e 

2m* 



(4) 



where denotes the lower subband edge and to* is the 
effective mass of the electron averaged over the extension 
of the wave function in well and barrier. 
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FIG. 1: Mixed momentum and real-space picture of a two- 
level-system that serves as a simple model for a diagonal 
intersubband transition. Tunneling into a virtual interme- 
diate state (dotted) as well as photon-assisted tunneling is 
expressed by a transfer matrix element hflij. Relaxation is 
assumed to take place within each subband only. 



A. Coherences 

Using the equation of motion of the density matrix 
ihdtp — [HtP] and separating the diagonal and non- 
diagonal part with respect to the parallel momentum k, k' 
according to p k \, = S kk > p\ 3 + (1 - S k k')p kk , , we obtain 
two coupled equations with four terms each, which deter- 
mine the time evolution of the system. Since the coherent 
term H is diagonal with respect to the in-plane momen- 
tum and the scattering term V is purely non-diagonal, 
the commutators that determine the diagonal and non- 
diagonal part of the density matrix are evaluated as 



ifid t p k 3 



pTK t ) 



(5) 



Thus, electrons are allowed to tunnel between the sub- 
band state i and j by means of the momentum conserving 
matrix elements 7if2y , in each of which they are possibly 
scattered out of a virtual intermediate state by an intra- 
well relaxation process V kk , as depicted in Fig. |l|. For 
simplicity we restrict ourselves to elastic (impurity, e.g.) 
scattering within each subband and assume a parabolic 



ptk'H^ j ), 



(G) 



where the commutator of the scattering potential with 
the non-diagonal part of p has been neglected in the sec- 
ond equation (Born approximation). The steady-state 
values of the coherences of the density matrix / lJ , which 
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determine the transitions \ik) — > \jk), the current a: 
the absorption, are obtained from a Laplace average' 15 ' 
defined by 



s I dte~ st p(t) 



(7) 



and performing the Laplace limit s — > using the relation 
lim s ^o + (u> — is) -1 = P(l/u>) + iir5(u;) at the appropriate 
stage of the calculation. In this approach the populations 
f l k l of the density matrix are not accessible and appear 
in the resulting expression as external quantities. The 
Laplace average gives 



ifisfl 3 



rn 

+ Y,<y&f2i-f&v$) ( 8 ) 

m 
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In a first step, the non-diagonal part in equation (|9|) is 
approximated (cf. appendix A for details). Specifying 
to the assumptions of a two-level-system with intra-well 
scattering only and neglecting terms of higher order in 
corresponding to multiple tunneling processes, gives 
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which has to be placed into equation (g) for the diagonal 
part. Simplifying for intra-well scattering here and taking 
the Laplace limit yields (cf. appendix B) 
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where (/) denotes the approximated expression for the 
non-diagonal part in equation ([To|). Performing an en- 
semble average, i.e. dropping terms related to correlation 
effects in the scattering potential, we obtain 



(e ik - tjk)fk ~ ™f h 3 x 
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which agrees with the previous resultli. In contrast to 
the original treatment, we continue by neither neglect- 
ing the difference of the arguments in the ^-functions on 



the LHS nor omitting the second term on the RHS. The 
coherence associated with the transition \2k) — ► \lk) is 
obtained from 



transition broadc 



rlation difference 
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where we have used abbreviations for the scattering in- 
duced broadening of the transition 7J, = J2 k , 6(ei k > — 
e jk)\Vkk> I 2 ' the subband separation e = e 2k — fife 
and the in-plane momentum of the final state q± = 
hT 1 yj2m*(e k ± e) . The first term on the RHS, which 
contains the difference of populations between the two 
states, corresponds to the central result of Referenced, 
cited often as the original proposal of the quantum cas- 
cade laser. The second term, which has been discarded 
so far, contains differences in population within a sub- 
band. It is this term which will be responsible for the 
second order type of gain, leading to the characteristic 
negative differential conductivity and the dispersive gain 
profile in a superlattice, and a modified spectral shape of 
the gain in a quantum cascade laser. 



B. Current density 

The current density between two states spatially sep- 
arated by d — z 22 — zn is calculated from j = eTr(vf), 
where v — i/h[H, z] is the velocity operator and z is the 
position operator. The current is induced by the non- 
diagonal matrix elements of the velocity operator uy, 
which are given by Vij — iQij^Zjj — Zy) 4- eZij/K. By 
choice of the basis set, the contribution_af the dipole Zij 
is small compared to the tunneling termES and we obtain 



(14) 



Using the previous equation for the coherences f^ 2 and 
/I 1 and the current yields 

. _ ed\nn 21 \ 2 ^ i\(f 2 k 2 - fjlHjlifl fl l ) 
3 n ^ e 2 + W + 7, 2 ) 2 - [ } 

The current results from differences in population. In the 
following section, and in contrast to the original work, the 
differences are evaluated for non-equivalent /c-states in 
the respective subbands. To obtain the result of Kazari- 
nov and Suris, q± is set equal to k and a constant broad- 
ening 7 is used 

. ed\nn 21 \ 2 7 

h e l + 7^ 

With this approximation, the current density is solely 
driven by the density of excess electrons in either state 
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An = J2k(fk 2 ~fl 1 )- For a superlattice, this approxima- 
tion does predict the resonant current peaks that occur 
whenever ground and excited states align, but fails to 
account for the current between equivalent states in the 
Wannier-Stark ladder. 



C. Absorption and gain 

Optical properties are deduced from the high- 
frequency response to an additionally applied ac field. 
In the case of a photon-assisted tunneling transition, the 
Hamiltonian is supplemented by 
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for a vector potential A = (c/iui)f UJ e~ % ' JJt with an ampli- 
tude fuj of the high-frequency field. Noting the similar 
structure of 8H and the non-diagonal part of H, the cor- 
rections to the coherences in linear response Sf 21 and 
8fl 2 are related to fj; 1 and f^ 2 by 
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edfu, ,21 , c,12 _ ed fu ,12 

Jk,e-hu ancl °Tk,U) — *^ Jk,e+huj> K 1 ' ) 



Kid 



which are evaluated at an energy eihto instead of e due to 
the time dependence of the ac field. This relation reflects 
the similarity of tunneling and photon- assisted tunneling 
processes in a diagonal transition. The photon-induced 
current becomes 



Sj(u) « edj^ ^2x8/1% 
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(18) 



The high-frequency conductivity is related to the current 
by ct(cj) = d(j + Sj(u))/df u , and directly linked to the 
absorption (c/. appendix C for details) according to 
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The in-plane momenta of the final states are denoted 
by k± = h~ 1 ^2m*(ek ± (e — fiuj)). The preceding ex- 
pression contains the two gain mechanisms as limiting 
cases of a more general intersubband gain profile with 
a simple physical interpretation. Before we discuss the 
quantum-mechanical paths involved, the expression may 
be generalized to an arbitrarily located pair of subband 
states. 

To account for vertical as well as diagonal transi- 
tions the basis of eigenstates of the biased system - the 
Wannier-Stark basis for a superlattice - is chosen. Due to 
the assumption of intra-well scattering only, the dark cur- 
rent vanishes as the tunneling matrix element hflij is in- 
corporated in the extended wave functions. The photon- 
induced current, however, is then mediated by the dipole 



matrix element between the two subband states, which 
can no longer be considered as small in this basis. The 
non-diagonal part of the velocity operator is given by 
Vij = iezij. Since the operator of the high frequency field 
5H = —e/cvA does not change its purely non-diagonal 
structure, inspection of the previous equations and re- 
placement of idflij by iezij/h naturally extends the equa- 
tion for the gain profile to any intersubband transition 
and permits to omit the rather arbitrary distinction be- 
tween a diagonal and vertical transition 
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As will be shown in the following, expression (EOh allows 
a simple explanation of the gain mechanism in a superlat- 
tice and in a quantum cascade structure. It is instructive 
to rewrite the differences in populations as 

iiur - Hi) = ii (/f (i - Hi) - a - )) ,(2i) 

s * s / 

emission 1 2) — ► 1) absorption] 1) — >|2) 

which directly translate into the paths depicted in Fig. 
The two processes above relate the states \2k) and \lk + 
by the emission or absorption of a (non-resonant) pho- 
ton, huj 7^ e for k ^ fc + , assisted by relaxation within 
the lower state via 7^, which ensures momentum trans- 
fer. The second difference in equation (20) is interpreted 
accordingly, where the relaxation takes place within the 
upper state. 




FIG. 2: Possible quantum-mechanical paths: for an incident 
photon with energy hui 7^ e, absorption or stimulated emission 
may occur due to a non-resonant absorption or emission into 
an intermediate and a subsequent relaxation into the final 
state. Energy and momentum are conserved in this second 
order process. 

If one assumes constant and equal in-plane scattering 
times t = h/~/, a Fermi distribution with the same tem- 
perature T in each subband and chemical potentials [i\ 
and fj,2, respectively, the gain is analytically expressed as 
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(22) 
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where 5 = e — hw characterizes the off-resonant nature of 
the photon transition and 9(5) = 9(5) + 6{-5)e- s/k » T 
reflects the asymmetry between "too small" and "too 
large" photons with regard to the resonant transition. 



III. RESULTS 

In this part the theoretical model is evaluated and in- 
terpreted with respect to the gain profile of a superlattice 
and quantum cascade laser. 



Bloch oscillator 



In the case of a superlattice, populations and scatter- 
ing times are equal for symmetry reasons, f^ 1 = f% 2 and 
7fc = Ik- The characteristic negative differential con- 
ductivity of the current- voltage characteristic j(F) is re- 
covered from equation (|l^) if one identifies the subband 
spacing with the field drop per period, e = eFd, where 
F is the applied electric field. The current density reads 



m = 



edA 2 
Afi 



E 



7fc(/fc- ~ fk+) 

(eFd) 2 + (2 7fc )2 



(23) 



where we omitted the state indices and introduced the 
miniband width via A s=s AMl\2- The populations may be 
described by thermal distributions, either by Fermi-Dirac 
or Boltzmann statistics. The current-voltage character- 
istic resembles the Esaki-Tsu characteristic aud agrees 
quantitatively with the result by-Wacker et al£3 within 
the sequential tunneling picturellll valid for weakly cou- 
pled superlattices. This assumption is implicit in the 
present approach, as we allow only for next-neighbor in- 
teraction and multiple tunneling is excluded, correspond- 
ing to a limited coherence of spatially extended states. 

Rewriting the gain profile of equation ( po| ) specifically 
for a superlattice yields 
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(24) 



Note, as the miniband width of a superlattice_A, and 
the dipole matrix element, z^j, are related byli3 — 
dA/AeFd in the Wannier-Stark basis, the diagonal ex- 
pression of equation (|l9|) and the general expression of 
equation ( |2Cj ) for a(uo) provide identical results. At reso- 
nance, the incoming photon provokes transitions between 
equivalent states, k± = k, and absorption and emission 
balance each other as expected from a system with no 
population inversion, An = X)fc(/ 22 — /fc 1 ) = 0- ^ n the 
case of photons with energy Tiuj < e the lower state in- 
volved in this second order transitions will be less oc- 
cupied than the upper state, leading to an asymmetry 
between emission and absorption in favor of gain. In con- 
trast, for a photon energy exceeding the subband spacing 
absorption occurs. As illustrated in Fig. |^ equation ( p4| ) 



recovers the dispersive shape of the Bloch gain - giv- 
ing rise to absorption above and (stimulated) emission 
below the field-dependent Bloch frequency ujf, — eFd/h. 
Note, we have neglected any particularity of the actual 
scattering processes here. A detailed investigation in the 
framework of second order perturbation theoryE3 reveals 
a complex interplay of population effects and the influ- 
ence of the momentum transfer for the relaxation pro- 
cesses in a superlattice. 



peak gain 



Wannier-Stark 
spacing s=eFd, 
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T=77K ; 




photon energy hco in meV 



FIG. 3: Bloch oscillator at e = eFd — 15meV: a dispersive 
gain contribution arises at An = from second order pro- 
cesses by non-resonant photon emission/absorption followed 
by scattering events that ensure conservation of momentum. 
Whereas stimulated emission is predicted for %u> < e (cf. In- 
set: path of stimulated emission), absorption dominates for 
fih) > e. A constant in-plane relaxation time ofrs 0.2 ps 
is used. Parameters of the GaAs/AlAs-superlatticeE3: A = 



20.3meV, d = 5.1nm, 



We compare the Bloch gain derived within the present 
approach with the results obtained from the standard 
model based on semi-classical calculations. In the semi- 
classical approach!] the Boltzmann equation is solved in 
the relaxation time approximation for the distribution 
function f(k z ,k\\) of miniband electrons subject to an 
external dc and ac field F(t) = F + F w cos(wi), where 
F u <C F. In the case of a Maxwell distribution!] this 
yields 
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in the single relaxation time approximation and 
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+ (1 — ibJT e )(\ — lUJTp) 



for the improved two relaxation time approximation 
given by Ignatov et al., where distinct momentum and 
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energy relaxation times t p , r e are wed and which agrees 
with detailed Monte-Carlo studiestil. The ratio of Bessel 
functions contains the temperature dependence for a non- 
degenerate electron gas. Fig. || shows a comparison 
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FIG. 4: Semi-classical (dotted line) vs quantum-mechanical 
results (full line) for the absorption in a superlattice for differ- 
ent temperatures T. We assume a temperature-independent 
scattering time r = 0.2 ps in the quantum-mechanical model 
and set Tk >e — r in the semi-classical model. In the 
semi-classical picture the peak gain scales with the ratio 
I 1 (A/2k B T)/I (A/2k B T). The quantum mechanical gain 
profile exhibits an additional narrowing with lower temper- 
ature. 

of the semi-classical results and the quantum mechan- 
ical predictions for the same constant relaxation time, 
r = h/j = 0.2 ps, at different temperatures T. No in- 
dependent parameters are used. The two approaches 
agree remarkably well at high temperatures in the semi- 
classical limit eFd < A. The narrowing of the Bloch gain 
profile with lower temperature, compared to the semi- 
classical curve, reflects an explicit influence of the elec- 
tron distribution within the subband. This influence is 
absent in the semi-classical treatment, regardless of the 
approximation for the distribution function. In real de- 
vices, however, the electron temperature reaches 100 K 
and above, if the superlattice is biased beyond the onset 
of Bloch oscillations, according to a self-consistent theo- 
retical analysis of the in=plane distribution function in the 
Wannier-Stark pictureE3. Still, at electron temperatures 
of about 77 K, the considered superlattice, e.g., with a 
sheet density of = 10 10 cm" 2 exhibits a peak mate- 
rial gain of about 15 cm -1 , which exceeds the estimated 
value for waveguide losses in the terahertz rangeo. 



Bloch type contribution add up as shown in Fig. ^. For 
a negligible lower state population, An/n « 1, where 
n = J2k(fk 2 + /fc 1 )' equation (|2(]) is dominated by res- 
onant photon emission due to the population inversion 
between equivalent /c-states. It resembles the Lorentzian 
shaped inversion gain profile, linearly depending on the 
population inversion An. On the other hand, in the limit- 
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FIG. 5: The population inversion, An, leads to nearly sym- 
metric gain profile centered at hui = e (cf. left hand side of 
Inset, dashed line). In addition, the dispersive Bloch gain con- 
tributes also here (cf. right hand side of Inset, dotted line). 
The generalized intersubband gain consists of both contribu- 
tions. Its spectral shape becomes more and more asymmetric 
with decreasing population inversion. The sample parame- 
ters correspond to the terahertz quantum cascade lasersOu at 
about hco = 18.7 meV. A constant relaxation time r = 0.5 ps 
and populations n.2 = 3 • 10 9 cm" 2 , n\ = 1 ■ 10 9 cm -2 are used, 
An/n « 0.5. 

ing case of equal populations, the Bloch type contribution 
results in a dispersive gain profile as in the superlattice. 
In between, there is a smooth transition of the (usual) in- 
tersubband gain profile to the dispersive Bloch gain with 
decreasing An/n as shown in Fig. ^. 

Thus, equation (|()|) for the gain profile states, that 
there is a dispersive contribution to the gain profile in 
any intersubband transition, with a rising significance for 
An/n tending to zero. This result implies two predic- 
tions for an intersubband emitter such as the quantum 
cascade laser. First, the gain does not linearly depend on 
An, but there is a non-negligible intersubband gain even 
without a population inversion, scaling approximately 
linearly with the electron density n in the system. Sec- 
ondly, above threshold, the peak gain and, thus, the laser 
signal is expected to shift to lower energies with increas- 
ing current or temperature, as the ratio An/n generally 
decreases. 



B. Quantum cascade laser 

In the quantum cascade laser, the populations of the 
respective subband state depend on the design as well as 
current and temperature. Then, the inversion gain and a 



IV. DISCUSSION 

In sections II and III of this paper, the homogeneous 
broadening 7^, of intrasubband relaxation processes has 
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FIG. 6: Evolution of the generalized gain profile from the in- 
version gain - linearly depending on the population inversion 
- to the red-shifted Bloch gain with decreasing An/n, respec- 
tively increasing the lower state population n\ , while keeping 
the upper state density ri2 constant. 

been introduced, without specifying it in detail. For 
the numerical evaluations it was taken as fc-independent 
quantity. As mentioned before, a realistic calculation 
for microscopic interaction processes in a superlattice 
will be given elsewhereEJ. Moreover, the effect of in- 
homegeneous level broadening due to interface rough- 
ness has not yet been discussed. Inhomogeneity can be 
considered in a simplified ajpneoach such as the "local 
quantum-mechanical model"c3'E3. This model assumes a 
"global quasi- Fermi- level" in each subband, independent 
of the in-plane position. Within this model the disper- 
sive shape of the Bloch gain is not obscured by inhomo- 
geneous broadening. This holds true even in a diagonal 
structureE-3, where the subband fluctuations are not cor- 
related, though the line width will be determined by in- 
homogeneous broadeningLJ. Furthermore, intersubband 
plasmons are known to alter the line shape of the inter- 
subband transition and to cause a blue-shift of the in- 
tersubbancLxesonance due to dynamical screening of the 
dipole fielcO. However, the subband states of a Bloch 
oscillator or quantum cascade laser are generally weakly 
populated w 10 10 cm~ 2 compared to the onset of the col- 
lective phenomenal 2 !! beyond some 10 n cm~ 2 . 

In conclusion, the proposed model provides a uni- 
fied description of optical transitions between two- 
dimensional subbands. The upper and lower subband 



can be either of the same kind (superlattice) or of differ- 
ent kind (quantum cascade laser structure). 

In a superlattice a population inversion between equiv- 
alent states differing in energy by e = eFd cannot occur, 
as J 22 — fl 1 , due to the translational symmetry of the 
system. Consequently, resonant stimulated photon emis- 
sion processes are exactly balanced by the corresponding 
absorption processes. However, non-resonant second or- 
der processes exhibit gain for fuj < e, whereas net ab- 
sorption occurs for hco > e. This inversionless gain at 
Tiuj < e represents the quantum mechanical analogue 
to the Bloch gain predicted by semi-classical models, 
which had not been described previously. The quantum- 
mechanical approach agrees remarkably well with the 
semi-classical results in the high-temperature limit. In 
contrast to the miniband pictures, it provides an easily 
conceivable interpretation of the gain mechanism. 

In a quantum cascade laser structure, upper and lower 
state generally exhibit a different population, and in the 
ideal case, an inverted population. If the inversion de- 
creases, the quasi-symmetric gain spectrum at a high de- 
gree of inversion, where An/n w I, evolves to the dis- 
persive Bloch gain, for An/n « and below. In contrast 
to the peak inversion gain, which does not depend on 
temperature and decreases with scattering, the latter de- 
creases with temperature as the differences in occupation 
between initial and final states diminish, but increases 
with more frequent scattering processes. 

The theory predicts amplification without inversion 
below the intersubband resonance of two broadened 
states. The peak gain in any amplified intersubband 
transition relying on a poor population inversion, i.e. 
An/n « 0, exhibits a red-shift of the order of the level 
broadening 7 with respect to the transition energy. This 
dispersive gain contribution, that is responsible for the 
Bloch oscillator effect and which escaped observation so 
far, is expected to be experimentally accessible in a quan- 
tum cascade structure by a search for the attributed red- 
shift. 
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Appendix: Mathematical Details 

A. Non-diagonal part ft?., 

Equation (^) governs the dynamics of the non-diagonal 
part in k,k' 

rn 

+ ^ H ^fw-mK j )- (27) 

Assuming intra- well scattering only, the second and third 
term on the RHS yield 

E(\Hm cmj fimj/mj\ ^ t/ m rij fij-irjj 
\ v kk' Jk> Ik v kk< ) ~ V kk' Ik' Ik V kk' 

m 

El Trim rmj f im Tjmj\ _ 
\ n k Ikk' ~ Ikk' U k' ) — 

n I 

= (e ik -e j „)fj?„+hil ij {f i k i,-f&,). 

Neglecting the non-diagonal matrix element p kk , (0) and 
taking the Laplace average equation ( p7|) gives 



Jkk' 



V- 



1 



€ik ^jk' 



i7r5(e tk - cjk') 



The non-diagonal f kk , still depends on f kk , — f kk , . The 
coherences between states k and k' within the same sub- 
band are derived from the special version for i = j of 
equation ( p7| ) 

m 

+ J2(Htfkt'-ftk'H^), (29) 

771 

where the second and third term are evaluated as 



E 

m 


(Trim rmi 
\ v kk' Jk' ~ 


E 


( Tjirn rmi 
\ u k Jkk' — 




~ \^ik ^ik' 







fjl 



and terms of higher order in the tunneling matrix ele- 
ment corresponding to multiple tunneling processes are 
neglected. Taking the Laplace average yields 



( 28 ) & 



-V- 



1 



(30) 
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+VK5{e ik -e ik ,)VMm-ft) 



( fii f ii \ 

v kk'\Jk Jk'l 



The last term vanishes as either the (5-function or the 
difference in populations is zero. Placing the approxima- 

I 



tions for f kk , and f kk , in equation ( f27| ) and neglecting 
the principal value yields 



fkk' = inSfek -e jk >) HVLi 



\ri3 I fjj _ fjj \ \rii I fii 
v kk'\Jk Jk'l v kk'\Jk 



e jk e jk' 



f-ik — f-ik' 



+v k i,f$-n>v h 



A/,-' 



r 



(31) 



B. Diagonal part f k 

Equation (^) determines the dynamics of the diagonal 
part in k 

ihsf k j = ihs$(0) + ^(Htf^ - fk m K ] ) 

rn 

m,k' 

where (/) denotes the previous approximations of the 
non-diagonal part. The second term on the RHS is given 
by 

E (Hrfr-fk n Hr) = 

rn 

= (tik-e jk )ft j +hn ij (f» -ft) 
Performing the Laplace limit s — > we obtain 

(eik-e jk )f k j =nn ij (fi i -f k ij ) (33) 

-Ew& (frL-m-vii) 



If one assumes no correlation between scattering events 
due to the impurities in different wells, i.e. were to drop 
terms containing the product V kk ,V kk , for i ^ j, the prod- 
uct of scattering potentials and the approximated non- 
diagonal part becomes 



V k %(f)l k = m5{e lk , - e ]k )\V k %{ 



Af k 



and similarly 



(f\ij yp 
\J >kk' v k'k 



i^{e ik -e jk/ )\V^ k ,\ 2 



£jk €jk' 



Rewriting equation ( |34| ) and sorting terms finally leads 
to an equation for the relevant coherences between the 
two states, that determine transport properties such as 
the current density 



transition broadening 



population inversion— >/« 



(e ik - e jk )f k j - i*fj> E ~ I' + " *ik>)\V& I 2 = I^.AJT f k j ) 



(34) 



m> 



iTrE^^ - ejkW^-^-ifP - /») + S(e tk - e jk ,)\V» \ 2 (f» - /«) 



€ik' ^ik 



€j k - €j k > 



Bloch type contribution^/ 6 
I 

It corresponds to the results of Kazarinov and SurisEl 
after performing a non-correlated ensemble average on 

their equation and specifying on a two-level system. The ac _ field induced coherence is given by the trans- 

formation 



C. Absorption 



" huJ Jk,e-huj 



(35) 
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which yields like the two contributions on the RHS of equation (|3j 

add up in respect to the current density. We obtain 

°Jk — 



(e-^)( e -^-z( 7 2 + 7 i)) J' 

As each Sf — S(f qc + f ho ) the absorption consists of two 
contributions 

a(Lu) = a qc {uj) + a bo {u}), (36) 



a qc (co) 



e 2 d 2 \n 21 \ 

£on r cut 



E 



ft**- 1 



\ 



(e + M 2 + (^ fc -1 ) 2 (e - M a + (^fc 1 ) 2 



V 



(/f - A 11 ) 



main contribution 



/ 



(37) 



and 



\ 



b e 2 rf 2 |» 21 | 2 ^ 
e n r ca; ^ 



itift 2 ~ ft 2 ) -ilUll- ft 1 ) ittfF- - ft 2 ) - itUtl - ft 1 ) 



(e-M^CftTfe 1 ) 2 



main contribution 



(38) 



where we set fc± = h 1 y // 2m*(ek ± (e — huj)), l± — 

h- 1 y /2m*(e k ± (e + M) and ftr^ 1 = 7fc + t|. The first 
term, a 9C (u;), depends on the difference in population of 
equivalent fc-states and yields the usual Lorentzian line 
shape for the gain profile in case of a population inver- 
sion. The second contribution, a bo (uj), contains differ- 
ences of populations between different fc-states within the 
respective subband state and will be discussed later. 

In the following, we will omit the non-resonant con- 
tribution. The contribution accounts for the situation 
where the upper state lies below the lower state. How- 
ever, it is not neglected in the numerical calculations as it 
is important to prevent the divergence at oj = 0, which al- 
ters the line-shape for huj ~ 0(7), i.e. in the far-infrared 
or terahertz regime. If one regards the resonant contri- 



bution only, the inversion gain reads 

qc( , = e 2 d 2 |»2i| 2 V Ivr^jfl-Jf) 
1 ' eorircu Y ( e - hco) 2 + (Ht^ 1 ) 2 

whereas the Bloch type contribution gives 

6 = e 2 rf 2 |» 21 | 2 v ll(ft 1 + -ft 1 )-lt(ft 2 --ft 2 ) 
{ ' e n r cij ^ ( e _^)2 + (^ r -i)2 

The difficulty in assigning a path of transitions of the 
electron to the latter expression for ot° (u>) is resolved by 
adding both contributions (c/. equation (fly)). 



